Inelastic neutron scattering study of phonons and magnetic 

excitations in LaCoC>3 

Y. Kobayashil*] Thant Sin Naing, M. Suzuki, M. Akimitsu, and K. Asai 

Department of Applied Physics and Chemistry, 
The University of Electro-Communications, Chofu, Tokyo 182-8585, Japan 

K. Yamada 

Institute for Materials Research, Tohoku University, Sendai, 980-8577, Japan 

J. Akimitsu 

Department of Physics and Mathematics, 
Aoyama Gakuin University, Sagamihara 229-8551, Japan 

P. Manuel 

ISIS Facility, CCLRC Rutherford Appleton Laboratory, Didcot OX110QX, UK 

J. M. Tranquada and G. Shirane 
Brookhaven National Laboratory, New York 11973 USA 
(Dated: February 2, 2008) 



1 



Abstract 

We have investigated the phonon and the magnetic excitations in LaCoOs by inelastic neutron 
scattering measurements. The acoustic phonon dispersions show some characteristic features of 
the folded Brillouin zone (BZ) for the rhomb ohedr ally distorted perovskite structure containing 
two chemical formula units of LaCoOs in the unit cell. We observed two transverse optical (TO) 
phonon branches along (<5 S S), consistent with previously reported Raman active E g modes which 
show remarkable softening associated with the spin-state transition [Ishikawa et al, (Phys. Rev. 
Lett. 93 (2004) 136401.)]. We found that the softening takes place in the TO mode over the 
whole BZ. In contrast, the acoustic phonons show no anomalous softening associated with the 
spin-state transition. The low-energy paramagnetic scattering at 8 K is weak, increasing towards 
a maximum at E > 15 meV, consistent with excitation of the nonmagnetic low-spin to magnetic 
intermediate-spin state of Co 3+ ions. 
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I. INTRODUCTION 



Lanthanum cobalt oxide LaCoC>3 exhibits magnetic-electronic transitions around 100 and 
500 K . 1 ! 2 ! 3 ! 4 ! 5 ! 6 In order to explain these two transitions consistently, a two-stage spin-state 
transition model has been proposed. Within the model, the 100 K transition is associated 
with the thermal excitation of Co 3+ ions from the low spin (LS; 5 = 0) ground state to an 
intermediate spin (IS; 5 = 1) state, whereas the 500 K transition corresponds to a crossover 
from the IS state to a mixed state of IS and the high spin (HS; 5 = 2) state.— Reducing 
the volume either by hydrostatic pressure^ or by subsitution of rare-earth elements,—^ with 
a smaller ionic radius, for La causes the spin-state transitions to move to higher temper- 
ature, and makes the two transitions less distiguishable. Korotin et a/.— were the first to 
show theoretically that the IS state is energetically close to the LS state due to a large 
hybridization between Co-3<i and 0-2p orbitals. In a purely ionic picture, however, the 
ground state of Co 3+ ions is either the LS state (A c f > A ex ) or the HS state (A c j < A ex ), 
where A c f and A ex are the crystal-field splitting and the intra-atomic exchange interaction, 
respectively . 2 '^ 1 ^ Recently, it has been argued that the IS state might be understood by 
including spin-orbit effects in the purely ionic model,™ and comparisons have been made 
with electron-spin-resonance measurements.— Thus, the nature of the IS-state Co 3+ is still 
controversial. 

The 100-K spin-state transition couples with phonons. In fact, lattice anomalies cor- 
related with the spin-state transition have been reported. The lattice volume shows an 
anomalous expansion of about 1% around 100 K.— The elastic constant measured by lon- 
gitudinal ultrasound shows a dip of about 30% in magnitude around 100 K, suggesting a 
softening associated with coexistence of the LS and the IS-states of Co 3+ .— In addition, it is 
expected that the hybridization-induced IS-state of the Co 3+ ion should have an associated 
Jahn- Teller (JT) distortion, since the e g orbital is partially occupied. Several experiments 
suggest the orbital ordering due to the JT distortion. Yamaguchi et alM have reported an 
anomalous temperature dependence of the intensity and the splitting of the phonon modes 
with the energies of 35, 50 and 73 meV based on an infrared spectroscopic study, and pro- 
posed a local JT distortion resulting from the appearance of the IS-state Co 3+ . Louca and 
Sarrao have proposed a local static JT distortion based on the pair-density-function anal- 
ysis of the neutron scattering in the paramagnetic insulating phase of Lai-zSr^CoOa.— A 
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structural analysis using synchrotron X-ray diffraction, performed by Maris et aZ.,— infers 
an alternation of short and long Co-0 bonds accompanied with a subsequent change from 
rhombohedral R3c to monoclinic 12/ a for 20-300 K, suggesting the presence of the e g orbital 
ordering on the IS-state Co 3+ with a cooperative JT distortion of the Q2 type.— The JT 
distortion increases sharply around 70 K with increasing temperature. Recently, Ishikawa 
et al. have reported an energy shift and broadening associated with the population of the 
IS-state Co 3+ ions for the phonons of the E g rotational mode of O atoms at 11 meV, the E g 
vibration mode of La atoms at 21 meV, and the A\ g rotational mode of O atoms at 32 meV 
(energies are the values at 5 K) based on the Raman scattering experiment.— In addition, 
they have reported the appearance of the satellite peaks associated with the IS-state at the 
energies of 19, 45 and 84 meV with increasing temperature. Based on the results, they 
supported the orbital ordering due to the cooperative JT distortion proposed by Maris et 
a/.— These studies show that the investigation of the phonons is of paramount importance 
in understanding the spin-state transition. 

It is important to clarify the phonon dispersion over the Brillouin zone in order to inves- 
tigate the coupling between the phonons and the spin-state transition. The previous studies 
suggesting the lattice anomal y 6ll3ll4ll5ll6i:L8 probe only the zone center (T-point) of the Bril- 
louin zone. Only neutron scattering experiments using a single crystalline sample affords 
information over the Brillouin zone. Louca et al. have reported the q-integrated phonon 
spectra by neutron scattering study but they showed no information on the q dependence 
of the phonons.^ The phonon dispersion in the rhombohedral phase of the perovskite oxide 
over the whole Brillouin zone has not been reported to our knowledge although many studies 
of the soft-phonon mode near the structural phase transition in the cubic phase of LaA103 
have been reported . 19 > 20 

In this paper, we report the inelastic neutron scattering study on LaCoOs single crystal in 
order to determine the phonon dispersion and find the phonon anomaly associated with the 
spin-state transition of LaCoOs. In addition to the phonon scattering, we have investigated 
the inelastic magnetic neutron scattering associated with the energy gap between the LS 
and the IS state, which has been reported to be about 16 meV . 12 ' 21 
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II. SAMPLE PREPARATION AND EXPERIMENTAL PROCEDURE 



Single crystals of LaCoC>3 were grown with a lamp-image floating zone furnace by melting 
polycrystalline samples prepared by a solid-state- reaction of predried La2C>3 and CoO.— 
Laue photographs confirmed that each sample is a proper single crystal. LaCo03 has a 
rhombohedrally distorted perovskite structure with the crystal space group of R3c which 
contains two chemical formula units. Since the crystal contains four twins, each of which 
has the principal axis parallel to one of the (lll)-axes of the pseudo-cubic unit cell, and 
the rhombohedral distortion of the crystal from the cubic structure is small, we use the 
crystallographic indices based on the pseudo-cubic cell (a = 3.83 A and a = 90.6° at 295 
K— ) containing one molecule of LaCoC>3. We must mention that the R-point of the pseudo- 
cubic Brillouin zone (0.5 0.5 0.5) is the zone center (r-point) of the rhombohedral one due 
to the Brillouin zone folding. A rod of the single crystal with about 5 mm in diameter and 
about 40 mm in length and an assembly of similar four rods were used for the measurement 
at Brookhaven National Laboratory (BNL) and JAERI, respectively. The crystals were 
oriented with their [011] directions vertical. The misalignment of the four crystals with 
respect to each other is within 1°. 

The inelastic neutron scattering measurements on a triple-axis spectrometer (TAS) H7 
at BNL were performed in 1996. After preliminary measurements on PRISMA at ISIS in 
Rutherford Appleton Laboratory in 2003, the detailed inelastic neutron scattering experi- 
ments were performed on a triple-axis spectrometer TOPAN at JRR-3M in JAERI, Tokai. 
Both TAS experiments were performed with fixed final neutron energies (Ef) of 30.5 and 
14.7 meV, but the horizontal collimation was 40'-40'-40'-80' for the BNL measurement and 
30'-100'-60'-blank and blank-60'-60'-60', respectively for the JAERI measurement. For the 
phonon scattering, the scans were performed to measure both the transverse and the longi- 
tudinal phonons with the propagating vectors (5, 5, 6), (5, 0, 0) and (0, 6, 5). For the magnetic 
excitation study, we have performed inelastic neutron scattering measurements around (1 
0) and around an equivalent reciprocal lattice point (3 0). 
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FIG. 1: The inelastic neutron scattering intensities per 1000 k monitor for LaCoOs at Q = 
(1 — 5, 1 + S, 1 + 5) at T = 8 and 200 K with Ef = 14.7 meV, after being corrected by the 
temperature factor. The solid lines are fits to the data using a Gaussian line shape. The arrows 
indicate the peak position of the acoustic phonons. The dotted lines are the background level. 

III. EXPERIMENTAL RESULTS 
A. Phonon scattering 

The inelastic neutron scattering intensities at Q = (1 — 5,1 + 5,1 + 5) are shown in 
Fig. ^ The data were taken with the fixed final neutron energy (Ef) of 14.7 meV at 
temperatures (T) of 8 and 200 K. The intensities have been corrected by the temperature 
factor (n + 1) = 1/[1 — exp(— E/ksT)]. Relatively sharp peaks are observed at about 6 
meV for 5 = 0.1 and about 10 meV for 5 = 0.2 and 0.3. We identify these phonons as 
the transverse acoustic (TA) phonons propagating along q = (5, 5, 5). We confirmed these 
phonons also at equivalent reciprocal lattice points Q = (1 + 5, 1 — 5, 1 — 5). We find that 
the energy shift for the acoustic phonons between 8 and 200 K is small for all 5 although the 
energy width at 200 K is broader than that at 8 K. The broad peaks observed in the higher 
energy region (E ~ 15 meV) for 5 = 0.2 and 0.3 are due to the optical phonon scattering 
described later. We cannot identify the small peak observed at around 12 meV for 5 = 0.2. 
We performed a similar measurement along Q = (1 + 5, 1 + 5, 1 + 5) and observed the 
longitudinal acoustic (LA) phonons propagating along q = (5, 5, 5) (not shown). 

The inelastic neutron scattering intensities for the transverse optical phonons at T = 8 
and 200 K are shown in Fig. El The data were taken at Q = (3 — 5,1 + 5,1 + 5) with 
Ef = 30.5 meV. At 8K at Q — (2.5,1.5,1.5) (5 = 0.5), two broad peaks were observed 
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FIG. 2: The inelastic neutron scattering intensities per 1200 k monitor at Q = (3 — 8, 1 + 8, 1 + 8) 
at T = 8 and 200 K with Ef = 30.5 meV, after being corrected by the temperature factor. The 
solid lines are fits to the data using a Gaussian line shape. The arrows indicate the peak position 
of the optical phonons. The dotted lines are the background level. 

at around E = 13 and 22 meV. We confirmed that the two peaks are due to the phonon 
scattering by comparing the scattering intensities measured at other equivalent reciprocal 
lattice points [Q = (1.5,0.5,0.5) and (2.5, 1.5, 1.5)]. We consider that these two phonons 
at the R-points for the pseudo-cubic cell are identical to the phonons of the E g rotational 
mode of O atoms (11 meV at 5 K) and E g vibration mode of La atoms (21 meV at 5 K) 
observed in a Raman scattering experiment by Ishikawa et al. 18 although the latter phonon 
is slightly different in energy. It should be noted that the Raman active modes for the 
rhombohedral R3c symmetry are at the R-point of the pseudo-cubic Brillouin zone since the 
atomic displacements of the modes are anti-phase for the nearest neighbor pairs along the 
[lll]-a3dsWW 

The optical phonons at 5 = 0.5 show a remarkable softening and broadening with in- 
creasing temperature from 8 K to 200 K. We measured the phonon spectra also at 120 K 
and found that the optical phonon spectra are almost identical to those at 200 K. The tem- 
perature dependence of the phonon energy is in agreement with the result of the Raman 
spectroscopy (see Fig. |3J). The broadening of the phonon peak of the O-rotational mode 
with increasing temperature is also consistent with the result of Raman scattering^. 

The phonon dispersion along q = (5, 5, 5) obtained in the present experiment is sum- 
marized in Fig. 0] (a). In addition to the LA and the TA branches, we see two transverse 
optical branches; the lower one continuing to the E g rotational mode of O atoms and the 
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FIG. 3: Temperature dependence of the phonon energy at Q = (2.50,1.50,1.50). The Raman 
data are from Ref. Q. The dashed lines are guides to eyes. 
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FIG. 4: The phonon dispersion curves of LaCoOs along (a) (6 5 5), (b) (0 5 5) and (c) (5 
0). The lines are guides to eyes. Inset in Fig. 0] (c) shows the trajectories for the scans. The 
hexagon and the rectangle are the (011) c section of the Brillouin zones for the rhombohedral and 
the pseudo-cubic unit cells, respectively. 

other to the E g vibration mode of La atoms at 5 = 0.5 . 18 i 23 i 24 We name the former and the 
latter branches TOl and T02, respectively. The phonon energies of both TOl and T02 
remarkably decrease with increasing temperature from 8 K to 200 K. The softening is most 
pronounced at 5 = 0.5, but it is not limited to that point. In contrast to the optical phonons, 
the softening in both LA and TA branches is negligible in all of the reciprocal lattice space 
investigated. 

We measured the phonon scattering at Q = (5 2 2) and (4+5 0) with Ef = 30.5 meV 
at 120 K to investigate the phonon dispersion along (5 0) for the TA and the LA modes, 
and Q = (2 5 5) with Ef = 14.7 meV at 250 K for the TA phonon mode along (0 5 5). The 
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phonon dispersion curves are shown in Figs. |U(b) and (c), respectively. The trajectories of 
the q in the reduced Brillouin zone for the rhombohedral crystal structure (R3c) containing 
two chemical formula units of LaCoC>3 are shown in the inset of Fig. 0] (c). 

B. Magnetic scattering 

The energy-dependent spectra at Qi = (1.04, 0.04, 0.04) and an equivalent reciprocal 
lattice point Q 2 = (2.96, 0.04, 0.04) measured at temperatures 8 and 120 K are shown 
in Fig. El We expect the signal at Qi to be dominated by magnetic scattering and that 
at Q2 to be dominated by phonons. Let us consider first the Q 2 data, in Fig. Efb) and 
(d). To fit these spectra, we assume that the scattering can be decomposed into three 
components: 1) gaussian tail from elastic scattering, 2) the two gaussian phonon peaks, 
and 3) a temperature- and energy-independent background. The magnetic scattering at Q 2 
should be small due to the magnetic form factor. In fitting the data, the background is 
constrained by the assumption that the temperature dependence of the phonon scattering 
is due to entirely to the Bose factor. The various fitted components are indicated in the 
figure. 

To extract the magnetic scattering from the spectra taken at Qi, we assume that the 
elastic tail and the background are the same as those determined at Q 2 . The phonon 
intensity is estimated by assuming that it varies as \Q\ 2 - The components so obtained are 
indicated in Fig.Efa) and (c). We assume that the difference between the measured data and 
the sum of the three components (tail, background, and phonons) is magnetic scattering. 
The deduced magnetic scattering at 8 K is shown in Fig. |U1 The significance of this signal, 
and the increased scattering at 120 K, will be discussed in the next section. 

IV. DISCUSSION 

A. Phonon scattering 

The features of the acoustic phonon dispersion along the principal axes of the reciprocal 
lattice space for the pseudo cubic unit cell are as follows. (1) For q = (5 5 5), both LA and 
TA phonon branches are folded about 5 = 0.25. (2) For q = (0 5 5), the energy of the TA 
branch increases with increasing 5 for 5 < 0.3, but levels off with E ~ 10 meV for larger 5. 

9 
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FIG. 5: The energy spectra of the inelastic neutron scattering intensities at Qi = (1-04, 0.04, 
0.04) (a) at 8 and (b) 120 K, and Q 2 = (2.96, 0.04, 0.04) (c) at 8 and (d) 120 K. The solid lines 
are the Bragg (or some incoherent) scattering tail. The dotted lines are the phonon component 
including the background intensity shown by the dash-dotted lines. 



f ig_6 . eps 



FIG. 6: The magnetic scattering intensity at Qi = (1.04, 0.04, 0.04) at 8 K. 

(3) For q = (5 0), the energies of the LA and TA branches increase monotonically with 
increasing 5. The features (1) and (2) can be understood in terms of the reduced Brillouin 
zone for the rhombohedral crystal structure of LaCoOa (see the inset of Fig. |H (c)). For q 
= (5 5 5), q crosses the zone boundary vertically at 5 = 0.25, and the dispersion should be 
symmetric about this point. For q = (0 S S), q arrives at the zone boundary with 5 = 0.375, 
and travels on the zone surface up to 5 = 0.5. The result shows that the phonon energy on 
the trajectory within the surface is not significant. The dispersions along (5 5 5) and (0 5 
8) are very different from those in the cubic phase of SrTiC>3. 25 

It has been reported that the longitudinal sound velocity along [111] is more than 10 % 
smaller at 200 K than at 8 K.— However, the temperature difference between 8 and 200 K 
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of the phonon energy found in the present study is substantially smaller than that (at most 
5%) for both the LA and the TA modes along (5 5 5) in the whole Brillouin zone, which 
means that the softening found in the ultrasonic measurements must be limited to a narrow 
region close to 5 = 0. 

The optical phonons show remarkable softening with increasing temperature in accor- 
dance with the Raman spectroscopy results.— A new finding of the present study is that 
the softening occurs not only at (0.5, 0.5, 0.5) but in the whole Brillouin zone along (5 5 
5) although it is most pronounced at 5 = 0.5. It should be noted that the E g rotational 
mode of O atoms (TOl at 5 = 0.5) and the E g vibration mode of La atoms (T02 at 5 
= 0.5) show substantial effects even though they are not especially sensitive to the force 
constant between Co and O atoms. More direct information about possible JT distortions 
is contained in phonons such as the quadrupolar mode found at much higher energies than 
those studied here. It would be quite interesting to probe that regime in a future study. 



B. Magnetic scattering 

The magnetic scattering at Q\ = (1.04, 0.04, 0.04) shown in Fig. |U] is weak at low 
energy and suggests a broad peak at ~ 20 meV. This energy scale is comparable to the 
energy gap between the LS and the IS states deduced from magnetization- (18 meV), lattice 
volume expansion^ (21 meV), NMR— (15.5 meV), and ESR— (12 meV) measurements. 
The similarity of these energies suggests that the increment of the magnetic scattering is 
associated with the excitation from the LS to the IS state of Co 3+ ions. In other words, the 
energy gap between the LS and the IS states are directly observed by the present inelastic 
neutron scattering. However, we regard this result as tentative because the energy region 
about ~ 15 meV is inevitably contaminated by the optical phonons, and furthermore we 
cannot obtain reliable temperature dependence of the magnitude of the excitation due to the 
increasing paramagnetic scattering arising from the fluctuation of the IS magnetic moments 
with increasing temperature. 

The difference of the neutron scattering between 8 and 120 K at Qi = (1.04, 0.04, 0.04) is 
shown in Fig. 0(a) along with the paramagnetic scattering intensity at (1.07, 0, 0) at 150 K 
measured previously with polarized neutrons. 3 Solid line represents the magnetic scattering 



function 



S(Q,u) oc x(Q)^r/(cj 2 + r 2 )]{l/[l-exp(-t,;/kT)]} with T = 4 meV, convolved 
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FIG. 7: The magnetic scattering intensity obtained by 1(120 K) - 1(8 K) measured at (a) Q\ and 
(b) Q2- Polarized neutron scattering data£ are also shown. The solid lines denote the resolution- 
convoluted scattering function^ The dotted line shows the phonon scattering intensity. 

with the instrumental resolution (r res = 5 meV). We see that the increment of the scattering 
can be ascribed to the fluctuation of the IS magnetic moments and is represented with the 
simple scattering function up to energies as high as 25 meV, although the solid line is only 
a guide to the eye. Here the increment of the phonon intensity between 8 and 120 K (shown 
by dotted line) is negligible. The difference of the neutron scattering between 8 and 120 K 
at Q 2 = (2.96, 0.04, 0.04) [Fig. 0(b)] is also interpreted by the paramagnetic scattering and 
increased intensity of the phonon scattering. Here, the former and the latter are denoted 
by the solid and dotted lines, respectively. The paramagnetic scattering intensity at Q2 is 
consistent with that at Qi after being normalized by the magnetic form factor. 

V. CONCLUSION 

We have investigated the phonons in LaCoOa using inelastic neutron scattering. The 
phonon dispersion curves for high-symmetry directions and E < 25 meV have been clarified. 
For (S S S), both LA and TA phonon branches show a symmetric shape about 5 = 0.25. For (0 
S S) , the phonon energy of the TA branch levels off for 5 0.3. These behaviors are consistent 
with the folded Brillouin zone for the rhombohedral perovskite structure containing two 
chemical formula units. The transverse optical branches continuing to the E g rotational 
mode of O atoms and the E g vibrational mode of La atoms show remarkable softening 
associated with the spin-state transition. The softening occurs over much of the Brillouin 
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zone. In contrast, the acoustic phonons show a much smaller change with temperature. The 
magnetic scattering intensity at around (1 0) increases to a maximum for E ^ 15 meV, 
ascribable to the excitation of Co 3+ from the LS to the IS spin-state. At higher temperatures, 
the paramagnetic scattering associated with the population of the IS-state Co 3+ ions was 
observed up to 25 meV. 
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